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O ' Abstract 
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The sum of all the ladder and rainbow diagrams in cj) theory near 6 dimen- 

vci . 

Q\ , sions leads to self-consistent higher order differential equations in coordinate 

^. 

space which are not particularly simple for arbitrary dimension D. We have 
D , now succeeded in solving these equations, expressing the results in terms of 

generalized hypergeometric functions; the expansion and representation of 
these functions can then be used to prove the absence of renormalization fac- 
tors which are transcendental for this theory and this topology to all orders in 
perturbation theory. The correct anomalous scaling dimensions of the Green 

functions are also obtained in the six-dimensional limit. 
ll.10.Gh, ll.10.Jj, ll.10.Kk 
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I. INTRODUCTION 



In a recent paper M we managed to derive closed forms for ladder corrections to self- 
energy graphs (rainbows) and vertices, in the context of dimensional renormalization. We 
only succeeded in carrying out this program for Yukawa couplings near four dimensions, 
although we did obtain the differential equations pertaining to the <ft 3 theory as well; but 
we were not able to solve the latter in simple terms. We have now managed to obtain 
closed expressions for 3 theory as well and wish to report the results here. The answers 
are indeed non-trivial and take the form of qF 3 functions, which perhaps explains why they 
had eluded us so far. Interestingly, the closed form results for Yukawa-type models lead to 
Bessel functions with curious indices and arguments; but as these can also be written as o-^i 
functions, the analogy with <p z is close after all. 

Given the exact form of the results, both for rainbows and ladders, we are able to test 
out Kreimer's |2| hypothesis about the connection between knot theory and renormalization 
theory with confidence, fully verifying that the renormalization factors for such topologies 
are indeed non-transcendental. At the same time we are able to determine the Z-factors 
to any given order in perturbation theory and show that in the D — > 6 limit, the correct 
anomalous dimensions of the Green functions do emerge, which is rather satisfying. The 
various Z-factors come out as poles in l/(D — 6) when the Green functions are expanded 
in the normal way as powers of the coupling constant, but the complete result produces the 
renormalized Green function to all orders in coupling for any dimension D. 

In the next section we treat the vertex diagrams, converting the differential equation 
for ladders into hypergeometric form. Upon picking the correct solution we are able to 
do two things: (i) establish that in the D — > 6 limit one arrives at the correct anomalous 
scaling factor for the vertex function, and (ii) obtain the Z-factors through a perturbative 
expansion of argument of the hypergeometric function. The case (i) is a bit tricky; it requires 
an asymptotic analysis, because the indices of the hypergeometric function as well as the 
argument diverge in the six-dimensional limit. The next section contains the analysis of the 
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rainbow graphs; the equations are similar to the vertex case, but different solutions must 
be selected, resulting in a different anomalous dimension. It is nevertheless true that the 
self-energy renormalization constant remains non-transcendental. A brief concluding section 
ends the paper. 



II. LADDER VERTEX DIAGRAMS 

We will only treat the massless case, since this is sufficient to specify the Z-factors once 
an external momentum scale is introduced. To further simplify the problem we shall consider 
the case where the vertex is at zero-momentum transfer, leaving just one external momentum 
p. The equation for the 1-particle irreducible vertex T, in the ladder approximation, thereby 
reduces to 

r,/ n „ d D q/(2n) D 

Letting T(p) = p 4 G(p), the equation can be Fourier-transformed into the coordinate-space 
equation for G, 

[d i -ig 2 A e (x)]G(x) = Z6 D (x), (2) 

where A c is the causal Feynman propagator for arbitrary dimension D. Since the coupling 
g is dimensionful when D ^ 6, it is convenient to introduce a mass scale \i and define a 
dimensionless coupling parameter a via, 

g 2 T(D/2 - 1) _ 4a(/2r) 6 - D 
Then, rotating to Euclidean space (r 2 = — x 2 ), the ladder vertex equation simplifies to 



d 2 D - 1 d 2 Aa(nr) 



G(r) = Z5 D (r). (3) 



dr 2 r dr r 4 

This is trivial to solve when D — 6, since it becomes homogeneous for r > and the 
appropriate solution is 
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G(r) oc r 6 ; b = -1 - \/§-2y/ATa, 



reducing to G(r) oc r~ 2 or r(p) = 1 in the free field case (a = 0); it represents a useful limit 
when analysing the full equation (3), to which we now turn. 

Let us define the scaling operator r = rj^. This allows us to rewrite the square of the 



d'Alembertian operator as 



d 1 



d 2 D - 1 d 



dr' A 



+ 



r dr 



r" 4 (0 r - 2)0 r (6 r + D- 4)(6 r + D-2). 



(4) 



Hence for r > the original equation (3) reduces to the simpler form, 



e p (e p — 2) (O p + D — 4)(© p + D-2)- Aap 6 ~ D G = 0, 



(5) 



where p = fir and 6 P is the corresponding scaling operator. Next, rescaling the argument 
to t — \.av^p~ x l v ', with v = 1/(D — 6), we obtain the hypergeometric equation: 



[6 t (e t + 2u)(Q t - 1 - 2v)(Q t - 1 - 4i/) - t] G = 0. 



(6) 



Being of fourth order, there are four linearly independent solutions 



0^3(61,62,63;*), 



t 1 - bl F 3 (2 - 61, 6 2 - 61 + 1, 63 - 61 + 1; t) 
t 1 - 62 i ? 3(2 - 6 2 , 63 - 6 2 + 1, 61 - 6 2 + 1; t), 

t 1 - fe oi ? 3 (2 - 63, 61 - 63 + 1, 62 - 63 + 1; f), 

where b± = 1 + 2z/, 6 2 = — 4z/, 63 = — 2z/. The appropriate solution, which near t — 
behaves as r 4_D when a = , is the last choice, namely 

G oc t 1+2l/ i 71 3 (2 + 2v, 2 + 4z/, 1 - 2v\ t). 
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Near r = this behaves like r~ 2 ~ l / v . Finally, renormalizing the Green function G to equal 
fi D ~ A when r = the scale we introduced previously for the coupling constant, and 

restoring the original variables, we end up with the exact result 

P In 2_ 9 4_ 1 j 2_. 4a(^r) 6 - p x 

- 31 6 - D ' e-D' 1-1- 6-£>' (6-g)* ) , . 

^ ' ~ T? (0 — 2 9 _ 4 i T 2 . 4a \ ■ V ' / 

0^3^ 6-D' Z 6^' 1 + 6-D'(6-Or 

To check that the poles in (D — 6) cancel out at any given order in perturbation theory, 
one simply expands the numerator and denominator in (7) to any particular power in the 
dimensionless coupling a and take the limit as D — > 6. For instance, to order a 3 , with a 
little work one arrives at, 

o, , , a 2 , , , . , . a 3 



r 2 G(r) -> 1 + - ln( fir) H \n(ur) (1 + 2 \n(ur)) H \n(ur) ( 9 + 6 ln(^r) + 4 ln 2 (/xr)) + 0(a 4 ). 

4 64 1536 v ' 

(8) 



It is most gratifying that this agrees perfectly with the expansion of the scaling index b 
obtained previously at D = 6. The most significant point is that there is no sign of a 
transcendental constant in the singularities of the perturbation expansion for 3 theory near 
6-dimensions, signifying that the renormalization constant Z is free of them, in agreement 
with the Kreimer hypothesis based on knot theory. 

One last (rather difficult) check on our work is to see what happens directly to (7) as D 
approaches 6, without having to invoke perturbation theory. For that an asymptotic analysis 
based on the method of steepest descent (see for example, de Bruijn [|J) is needed. We 
start by making use of the Barnes integral representation of the hypergeometric function, 

o a i, 2, 3, ) 2m J-ioo r(6i + z)T(b 2 + z)T(b 3 + z) V ; 
In our case the b arguments lead us to evaluate the integral 

KJ 2inJ-too r(2 + 4z/ + z)r(2 + 2z/ + ^)r(l -2v + z) 1 H J v; 



in the limit as v — > oo. We shall show that as a function of p, I u behaves like p 1 V 5 2 v / 3+a_ 
Remember that G(r) oc r- 2 ~ l / v I u (r). 



For the method of steepest descents, suppose we write I u as 

j MOO 

— / g v (z) exp[f„(z)] dz 

where g v is a "slowly varying" function. We see from (9) that all the poles of the integrand 
lie on the positive real axis. If ( is such that Re(£)< and f '„((,) = 0, then an approximate 
evaluation of I u is given by 



ag u {t)eMMC)]/p*\tt{Ql 

where 

a = expHaxg(/*(C))/2]. 

On applying the reflection formula for the gamma function || to both r(l — 2v) and T(l — 
2v + z) appearing in (9), we find that we can write 

sm{27ru) 

provided v is not an integer, and 

«p[M*)l = r(2 + 4,)r(2 + 2,)r(2,- g) r(-,) 

vuvK )\ r(2 + 4z/ + ^)r(2 + 2i/ + 2)r(2z/) v ; 

Since 

f v (z) = log(4az/ 4 ) - [V(2i/ -z)+ ij)(-z) + ^(2 + 2v + z)+ ^{2 + Au + z)], 

where ip denotes the psi (or digamma) function, we look for a zero at z = — £z/ say, where 
< £ < 2. Since we assume v ^> 1 and since for x ^> 1, ip(x) = logx + 0(l/x), we find that 
£ must satisfy the quartic 

£(£ + 2)(£-2)(£-4) =4a. 
The four solutions of this equation are 



£ = 1± v / 5±2 v / 4T 



and are all real if < a < 9/4. In particular we shall choose the zero j3 say in (0,2) which is 
closest to the origin; that is 

(3 = 1- \/5-2y4Ta. 

With this value of (3 we find 

fl(-(3v)~{l-(3){A + 2(3-(3 2 )/{av). 

Since in fact < f3 < 1, we have that arg/"(— (3u) = so that a — 1. Again, 

sm((2 + P)nu) 



and, after some algebra, 



g v {-(3v) = — — r p 

sin(27rz/) 



4V2^P(P + 2) l 16(2 + 0) , 2v 



eMM-P")\ = fl3/2l/1/2 [ ( 4 - m2 -/3) ] 6XP( ~ W ' 
approximately. Consequently, for v 1 but not an integer, we find 

F 3 (2+2,, 2+4,, 1-2,, 4avp ) ~ asin(27rz/) (1 _ ^ + 2/3 _ ^1/2 t (4 _ ^ _ ^ 
Using this asymptotic expansion, we obtain simply from eq.(7) that 

G(r) = ^-o-V^m^ (10 ) 

which is just the scaling behaviour at 6-dimensions which we were seeking. We have therefore 
fully verified the correctness of (7) in all the limits. The last step is to convert the answer 
to Minkowski space by making the familiar substitution r 2 — > — x 2 + ie. 



III. RAINBOW DIAGRAMS 

Let Ar(p) denote the renormalized propagator in rainbow approximation, so that 
p 2 A R (p) = 1 — T lR (p)/p 2 , where S R is the rainbow self-energy. The propagator obeys the 
integral equation in momentum space 
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p A R (p) = Zp + zg j {2n)D k2 , (11) 

where Z now refers to the wave- function renormalization constant. As always we convert 
this into an x-space differential equation, 

[«9 4 - ig 2 A c (x)]A R (x) = -Zd 2 5 D (x). (12) 

Interestingly, this is exactly the same equation as (2), apart from the right hand side, 
and it can therefore be converted into hypergeometric form by following the same steps as 
before. The only difference is that we should look for a different solution, because as g — > 0, 
A fl (p) - 1/p 2 , or A R (x) ~ (x 2 ) 1 - / 2 . 

A simple analysis shows the correct solution is 

t 1+4iy F 3 (2 + 4z/,2 + 6z/,l + 2u;t); t = Aav\^r)- l/u , v = 1/(D - 6), 

because this reduces to r 2 ~ D when a = 0. Actually we can solve (12) directly at D = 6 
when a ^ because it is a simple homogeneous equation leading to 

A R (x) oc 

and thereby determine the anomalous dimension from the exponent of r. Anyhow, the exact 
solution of the rainbow sum for any D and renormalized at r — 1/ fi is here obtained to be 

P (n 4_ o 6_ 1 2_. Mpr) 6 - D \ 

A M _ J-D ° f3lZ 6-D' Z 6-D' 1 6-D' (6-D) 4 n „n 

0-^31+ 6-D' Z 6-D' 1 6-D' (6-D) 4 / 

The numerator and denominator of (13), when expanded in powers of a will reproduce the 
(renormalized) perturbation series; to third order we find, in the limit as D — > 6, that all 
poles disappear and 

2 3 

r A A R (r) = l~ ln(/ir) + ^ ln(^r)(ll+6 ln(^r))- ln(/ir) (103+66 ln(//r)+12 In 2 (fir))+0 (a 4 

This coincides perfectly with the expansion of the scaling exponent at D — 6. 

Lastly we need to show that the D — > 6 limit of (13) collapses to the scaling behaviour 
found above, via an asymptotic analysis of the Barnes representation. We have indicated 



how this can be proven in the previous section and thus we skip the formal details to avoid 
boring the reader. The long and the short of the analysis is that no transcendentals enter 
into the above expressions for the self-energy (including their singularities, which are tied to 
the wave- function renormalization constant). These results confirm nicely the Kreimer |2| 
hypothesis that the Z-factors will be simple rationals for such topologies. 

IV. CONCLUSIONS 

We have succeeded in evaluating an all-orders solution of Green functions for ladder and 
rainbow diagrams for any dimension D in 3 theory; the results are non-trivial, involving 0-^3 
hypergeometric functions. We have demonstrated that, in the limit as D — > 6, the correct 
six- dimensional scaling behaviour (which can be separately worked out) is reproduced. One 
can likewise determine the exact solutions for massless bubble ladder exchange in 4 theory, 
because the equations are very similar: they are also of fourth order and can be converted 
into hypergeometric form too |7[]. 

More intriguing is the question of what happens when self-energy and ladder insertions 
are considered, so far as renormalization constants are concerned. A recent paper by Kreimer 
II has shown that such topologies with their disjoint divergences can produce transcendental 
Z in accordance with link diagrams that are of the (2,q) torus knot variety, where the highest 
q is determined by the loop number. It would be interesting to show this result without 
resorting to perturbation theory by summing all those graphs exactly, as we have done in this 
paper. (Kreimer cautions that multiplicative renormalization may screen his new findings.) 
The generalization to massive propagators || does not seem beyond the realms of possibility 
either, although it has a marginal bearing on Z-factors. 
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